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Using the black string between two branes as a model of a brane-world black hole, we compute 
the gravity wave perturbations and identify the features arising from the additional polarizations of 
the graviton. The standard four-dimensional gravitational wave signal acquires late-time oscillations 
due to massive modes of the graviton. The Fourier transform of these oscillations shows a series of 
spikes associated with the masses of the Kaluza-Klein modes, providing in principle a spectroscopic 
signature of extra dimensions. 
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Black holes are central to our understanding of gravity, 
and are expected to be key sources of gravity waves that 
should be detected by the current and upcoming genera- 
tion of experiments. Such a detection will not only con- 
firm the indirect evidence from binary pulsars for gravity 
waves, but will also also allow us to probe the properties 
of black holes and of gravity. In particular, this will open 
up a new window for testing modifications to general rel- 
ativity, such as those arising from quantum gravity the- 
ories. String theory for example predicts that spacetime 
has extra spatial dimensions, so that the gravitational 
field propagates in higher dimensions and has extra polar- 
izations. Recent developments in string theory indicate 
that Standard Model fields may be confined to a four- 
dimensional 'brane', while gravity propagates in the full 
'bulk' spacetime. This has spurred the development of 
brane-world models, such as Randall-Sundrum (RS) type 
models, which can be used to explore astrophysical pre- 
dictions 0] . RS type models have a five-dimensional bulk 
with negative cosmological constant, so that the metric 
is warped along the extra dimension. As a result, these 
models provide a new approach to the hierarchy problem, 
dimensional reduction and holography. 

The nature of black holes that form by gravitational 
collapse on an RS brane is only partly understood 0, 0| , 
and no exact solution is known for a black hole localized 
on one brane. If there is a second 'shadow' brane, the 
black string may be used to model large black holes on 
the visible brane, when the horizon on the brane is much 
greater than the extent of the horizon into the bulk Q • 
The black string reproduces the Schwarzschild metric on 
the visible brane but is not confined to the brane, since 
there is a line singularity at r = into the extra di- 
mension (see Fig. [y. The shadow brane can also in- 
troduce an infra-red cut-off to shut down the Gregory- 
Laflamme (GL) instability of the black string at long 
wavelengths If the shadow brane is close enough to 
the visible brane for a given black hole mass M, or if 
M is large enough for a given brane separation d, then 
GMe~ d / e /£ is above a positive critical value and the GL 
instability is removed (see below). This is the background 
model that we perturb. 
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FIG. 1: Schematic of the black string 

The brane separation is constrained from above by sta- 
bility requirements. It is also constrained from below. 
This follows since the brane separation is a massless de- 
gree of freedom, felt on the visible brane as a 'radion' 
field, so that the low-energy effective theory on the visi- 
ble brane is of Brans-Dicke type, with [f| Wbd = 3(e 2d / e — 
l)/2, where I is the bulk curvature radius. The shadow 
brane must be far enough away that its gravitational in- 
fluence on the visible brane is within observational lim- 
its. Solar system observations impose the lower limit 
w bd > 4 x 10 4 , so that d/t > 5. The allowed region 
in parameter space is shown in Fig. 2. Table-top tests of 
Newton's law impose the constraint I < 0.1 mm. This up- 
per limit defines a mass 0.1 mm/ 2 G ~ 10~ 7 M Q , so that 
for astrophysical black holes it follows that 2GM ^> I is 
easily satisfied: GM/£ > 10 7 (M/M Q ). 

The 5D black string is a solution of the Einstein equa- 
tions Gab = &£~ 2 9ab, with metric 

ds 2 = a 2 (- f dt 2 + f- 1 dr 2 + r 2 dtf) + dy 2 , (1) 

where f(r) = 1 - 2GM/r, a(y) = e^*, and the 
branes are at y — 0, d. Metric perturbations satisfy 

nh AB + X7 A V B h c c ~2V c \7 {A h B)c -8r 2 h AB = 0, (2) 

and gauge choices may be made to reduce the degrees 
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FIG. 2: The allowed region in parameter space. The stability 
boundary is well approximated by GM/£ = 0.1122e d// *. 

of freedom to 5, the number of polarizations of the 5D 
graviton. 

As in standard black hole perturbation theory 0, we 
use spherical harmonics Yj m and their gradients to con- 



struct a basis T 



(j,lm) 
AB 



for h,AB, with j — 1 ... 15. The 



metric perturbation is naturally split into two decoupled 
parts: polar [parity (—1)'] and axial [parity (— 
Suppressing the (Im) indices, the expansion coefficients 
are of the form C« = e luJT H^ (r, y), where r = t/GM 
and uj is a dimensionless frequency. We use two gauges, 
which are related by quadrature. We generalize the 4D 
Regge- Wheeler (RW) gauge by setting the coefficients of 
the most complicated harmonic tensors equal to zero, 
which makes the remaining coefficients gauge invariant. 
The RS gauge 0, given by h A A = = V B h A B = h Ay , 
has the advantage that the expansion coefficients are sep- 
arable: = e luJT W^\r)Z{y), where 



a 2 Z" - 2{aa')'Z = -m 2 Z. 



(3) 



Here, m is the effective mass on the visible brane of the 
Kaluza-Klein (KK) mode of the 5D graviton. The bulk 
wave functions Z(y) are the same as those for Minkowski 
branes (even though the bulk is no longer anti de Sitter), 
and the solution for m > is Z(y) = B^iml/a), where 
B2 is a linear combination of Bcsscl functions. 

Neglecting brane bending, the boundary conditions at 
y = , d in RS gauge are d y (a~ 2 h a p) = 0. When m = 0, 



the solution is Z 



For this zero-mode, the met- 



ric perturbations reduce to those of a 4D Schwarzschild 
metric, as expected. For m 0, the boundary conditions 
lead to a discrete tower of KK mass eigenvalues, 



m n = (z„,/£)e 



(4) 



where Yi(m„£) Ji(z„) = Ji(m n e)Y 1 (z n ) 0. The GL in- 
stability exists in the range < m < m cr it w 0A/GM 0. 
If mi > m cl jt, then the instability is avoided, and 
mi = m cr it defines the stability curve in Fig. [21 

Radial master equations. We generalize the standard 
4D analysis to find radial master equations for a reduced 



set of variables, for all classes of perturbations. In the 
5D case, the KK modes are governed by coupled mas- 
ter equations. These can be written as matrix- valued 
Schrodinger-like equations, 



d;r 



(5) 



where V is the potential matrix, x = p + 2 In (p/2 — I), 
and p — r/GM. For spherical (s-wave) polar perturba- 
tions, there is only one master variable, and 



V s _ + 6 - 18 /i 4 p e - 24 it V + 36 fi 2 p 3 

J~ p 3 (p 2 p 3 + 2) 2 ' 

(6) 

where p, = GMm. For I > 2 axial perturbations, there 
are two master variables and the potential matrix in the 
RW gauge is 



Va 2t 

T = tl1 ' 



l(l+l)p-6 p 2 
4p 3 1(1 + l)p 



(7) 



The diagonal elements of V a are identical to the 
RW potentials for spin-2 and spin-1 excitations of 4D 
Schwarzschild, plus a mass term. Massive KK modes 
of the 5D graviton include spin-1 perturbations. (For 
p = 0, the spin-1 contribution is pure gauge.) For the 
other types of perturbations, we find that the I = 1 and 
I > 2 polar perturbations are governed by two and three 
master variables respectively, while there is only one mas- 
ter variable for the axial p-wave. 

Stability. All master equations can be cast in the form 
C(p,uj)^ = 0, where £ is a linear differential operator. 
Stability holds if C is positive definite when the bound- 
ary conditions \f = arc imposed at x = ±00, and 
uj 2 < 0. We have verified the positivity of C an- 
alytically for all classes of axial perturbation. For the 
I > 1 polar case we performed a numerical search for so- 
lutions satisfying the above boundary conditions in the 
relevant region of (p,oj) parameter space. No solutions 
were found, suggesting that no instability exists. The op- 
posite is true for the s-wave case. The reason lies in the 
potential, which is shown in Fig. for different values of 
p. For small p, there is a potential well that can support 
a normalizable bound state with uj 2 < 0, which implies 
that L is not positive definite. This well vanishes for large 
p, suggesting that an instability exists for all modes with 
< p < pent- We used a standard WKB-inspired phase 
integral analysis to develop a new improved deriva- 
tion of the GL instability. Our approach leads to an ac- 
curate determination, p cr { t = 0.4301, in agreement with 
recent results from full numerical relativity . Stability 
is achieved if the mass of the first KK mode is such that 
GMmi > ^ crit (see Fig. 2). 

Gravity wave signals. Perturbations of a stable black 
string describe gravitational radiation that may in prin- 
ciple be observable by current or future detectors. The 
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total gravity wave signal at the observer (x — x b s ) is 
a superposition of the waveforms V'nl 7 ") associated with 
the mass eigenvalues m n of Eq. J3J). In Fig. 0] we 
present signals associated with the four lowest masses for 
a marginally stable black string, for I — 2 axial modes. 
These are obtained by re-introducing time dependence 
into the Schrodinger equation (|5j), via lo — > — i d / dr, and 
numerically integrating the resulting PDE with V = V a . 
We assume static Gaussian initial data centered about 
x = 50 and an observer at x \, s = 100. Because the mas- 
ter equations reduce to the 4D RW equation when /j = 0, 
the bottom curve exhibits the damped single-frequency 
(quasi-normal) ringing followed by a power-law tail famil- 
iar from the 4D analysis [jj. The massive mode signals 
are quite different |10| , showing much less damping than 
the massless mode, and late-time monochromatic oscil- 
lations instead of a featureless power-law tail. This is 
reminiscent of the behaviour of massive scalar fields in 
4D spherically symmetric spacetimes, where a detailed 
Green's function analysis shows 0| that the late-time 
signal is oc i~ 5 / 6 sin mt, irrespective of I. 

To reconstruct the total waveform 7/> to t for the I = 2 
axial measured by an observer on the visible 

brane, we need to specify an initial profile F for the per- 
turbation in the bulk. Using a conformal bulk coordinate 



and assuming that F and Z n are appropriately 



normalized, we have 



■0tot(V) = c nZ n {0)lp n {T) 



£FZ n d£, (8) 



where ^ n = 1. Hence, c 2 is the fractional energy in 
the n-mode. Massive KK modes introduce striking new 
features in the gravity wave signal that are potentially 
observable, as shown in Fig. [S] The strength of the new 
features is sensitive to the initial data F(£). Here we 
simply illustrate the possibilities via two forms of F(£) 
that correspond to very different situations. 

In Fig. the upper signal has d/i = 6 12], and the 
initial data F(£) corresponds to the zero-mode l/£ 2 with 
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FIG. 3: The s-wave potential for varying /j,. The critical po- 
tential is indicated by the heavy line. For fi < /U C rit, there is 
a bound state with u> 2 < 0, indicating an instability. 
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FIG. 4: Gravity wave signals associated with 3 massive KK 
modes and the 0-mode of I = 2 axial perturbations, for a 
marginally stable black string. The inset shows the zero-mode 
signal on a logarithmic scale. 

a cut-off applied one eighth of the distance to the shadow 
brane. This data is motivated by the results of numerical 
calculations of the bulk gravitational field around a small 
brane-localized black hole |2j. If the black hole mass is 
small enough to be in the unstable region of Fig. 2, we 
may think of the black hole as the endpoint of a GL in- 
stability 0, 0- The initial data is then thought of as 
arising from an encounter between this black hole and 
the much heavier black string. We find that in this case, 
the coefficients of the massive mode signals are of order 
10~ 4 , which suggests that the total waveform will only 
exhibit minor deviations from the GR prediction. This is 
certainly true for early times, where the signal is domi- 
nated by ordinary 4D Schwarzschild quasi-normal ringing 
(cf. the inset of Fig.QJ. However, as can be inferred from 
the individual massive waveforms in Fig. 4, the late-time 
signal is very different from the 4D case-it is both os- 
cillatory and very lightly damped. The inset shows the 
Fourier transform of the signal for late times. There are 
eight discrete peaks in the spectrum, corresponding to 
each of the eight non-zero mass modes. 

In the lower panel, we take d/l = 20 and Gaussian ini- 
tial data in the bulk centered halfway between the branes 
(if the center is moved closer to the shadow brane, the 
results are qualitatively similar). This could correspond 
to an event which 'mainly takes place in the bulk', such 
as the merger of two black strings. In this case, the zero 
mode waveform is dominated by the massive mode sig- 
nals, and the waveform is very different from the 4D case. 
The contributions from the odd modes m^, 777,5 and m*? 
are suppressed, which results in a late-time transform 
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FIG. 5: Composite I = 2 axial gravity wave signal (normal- 
ized) from the 9 lowest mass modes of a marginally stable 
black string. The upper panel corresponds to 'truncated zero- 
mode' bulk initial data with d/l = 6, the lower panel to Gaus- 
sian bulk initial data with d/l = 20. Insets show the Fourier 
transform of the respective waveforms for 300 < r < 1000. 
The percentage of energy in the 4D zero-mode part of the 
signal is 99.96% (upper) and 2.281 x 10~ 13 % (lower). 

with five principal peaks in the inset. These peaks in 
fact display a degree of fine structure, as shown in the 
blow-up, that deserves further investigation. 

Conclusions. Will the consideration of other types of 
perturbations or different multipoles in some way con- 
taminate the spectroscopic signal? We find that in all 
cases, as r — > oo, master variables with [i ^ behave as 
massive fields propagating on Schwarzschild spacetime, 
which is enough to guarantee the slowly-decaying oscil- 
lating tail. Massless modes for I > 2 are known to be 
subject to more damping than for I = 2, so if the mas- 
sive signal dominates the latter for late times, it will also 
dominate the former. 

Can the massive signal be resolved by realistic grav- 
ity wave detectors? There are two separate issues: the 
frequency of the massive modes and their relative am- 
plitude. From the formulas above, one finds that the 
discrete frequencies in the late-time tail are 

.^"-^'(Olmm/f) Hz. 



fn 



(9) 



For d/l > 5, z n is well approximated by J\{z n ) = 0. 
Note that unlike 4D black hole quasinormal modes, these 
frequencies are independent of the mass M. Taking I = 
0.1mm and d/l — 24, results in /i ~ 100 Hz, which is 
ideal for a LIGO detection of extra dimensions. For d/l — 
33, we find an optimal LISA signal, with fi ~ 0.01 Hz. 
For these parameters, the lower limits on black string 
masses are M > 100M Q and M > 10 6 M Q , respectively. 



In regard to the amplitudes, the most optimistic type 
of events involve bulk-based initial data, i.e., black string 
mergers. Since the zero-mode energy is small in these 
situations, we expect the strength of the massive mode 
oscillations to be comparable to the quasinormal ring- 
ing amplitude in the analogous 4D case. The situation 
is more delicate for brane based initial data, where the 
zero mode-tends to dominate the signal. We find that 
the relative massive amplitude tends to increase with de- 
creasing d/l. However, small d/l implies large /„, so the 
best prospect of seeing this type of signal lies with high 
frequency detectors. 

The discrete nature of the late-time Fourier transform 
of the black string waveform is the most important ob- 
servable feature of this model. Its detection would pro- 
vide clear evidence of extra dimensions. It would also give 
the direct spectroscopic measurement of the KK masses 
m„, which in turn provides information about d and I. 
Although we have analyzed a specific model of brane- 
world black holes, we expect that qualitatively similar 
features will arise for other models with compact extra 
dimensions, since they all have a discrete tower of mas- 
sive KK modes. Furthermore, since the massive modes 
travel below light-speed, there will be potentially observ- 
able time-delay in their arrival, which could be of or- 
der seconds or longer for distant sources. Also, the light 
damping of KK modes means that they may have a sig- 
nificant integrated contribution to the stochastic gravity 
wave background. 
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